HST Observations of the Wolf-Rayet Nebula NGC 6888 
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ABSTRACT 

We present Hubble Space Telescope WFPC2 images of a portion of the bright northeast rim 
of NGC 6888, the nebular shell physically associated with the Wolf-Rayet star HD 192163. The 
exposures are taken in the light of Ha A6563, [0 III] A5007, and [S II] AA6717, 6731. The images 
are used to constrain models of the ionization structure of nebular features. From these models we 
infer physical conditions within features, and estimate elemental abundances within the nebula. 
The results of our analysis, together with the degree of small scale inhomogeneity apparent in the 
images, call into question the assumptions underlying traditional methodologies for interpretation 
of nebular spectroscopy. 

The thermal pressure of photoionized clumps is higher than the inferred internal pressure of 
the shocked stellar wind, implying that the current physical conditions have changed significantly 
over less than a few thousand years. These results are discussed within the context of published 
three-wind evolutionary scenarios for the formation of the nebula. We also discuss the nature of 
a back-illuminated radiative shock driven into the cavity surrounding NGC 6888. 



Subject headings: H II regions - ISM: individual (NGC 
individual (HD 192163) - stars: mass loss 

1. Introduction 

Wolf-Rayet (WR) stars are thought to be the 
late stage of the evolution of stars more massive 
than ~ 30 M Q . WR stars are characterized as 
possessing significant stellar winds, large in both 
mass-loss rate and terminal velocity. Observations 
of 62 Galactic WR stars in the northern sky show 
that about half are associated with nebular emis- 
sion (Miller and Chu 1993), some of which are ring 
nebulae. Johnson and Hogg (1965) were the first 
to propose that these ring nebulae were formed by 
stellar wind interaction with the local interstellar 
medium. Attempts have been made to analyti- 
cally explain the formation and structure of such 
nebulae (Weaver et al. 1977; Van Buren 1986). 

One of the most prominent and best studied ex- 
amples of this class of object is NGC 6888. NGC 
6888 is associated with WR 136 (HD 192163), con- 
sidered to be a member of the Cygnus OBI associ- 
ation. Its large angular size (18' x 12') has made it 
an attractive object for investigation at a variety 



ISM: kinematics and dynamics - stars: 



of wavelengths. For a distance of 1450 pc (Wend- 
ker et al. 1975) the nebula has a physical size of 
7.6 pc x 5.0 pc. One of the first detailed stud- 
ies of NGC 6888 was that of Lozhinskaya (1970) 
who, assuming that the nebula is prolate ellipsoid, 
stated that the nebula is tilted ~ 30° out of the 
sky. Later studies (Parker 1978; Kwitter 1981) 
used emission-line ratios to show that the nebula is 
overabundant in nitrogen and helium and under- 
abundant in oxygen. Explanations for this enrich- 
ment include the transport of processed material 
from the stellar interior to the outer atmosphere 
prior to the onset of the RSG phase (Esteban and 
Vilchez 1992), and the mixing of interstellar and 
wind material (Kwitter 1981). The high level of 
nitrogen is consistent with the classification of WR 
136 as a WN6 star. An [O III] skin external to the 
Ha nebulosity has previously been observed (Du- 
four 1989; Mitra 1990), suggesting that the stellar 
wind and the expanding shell are interacting with 
the bubble outside the RSG shell. 

There have been a number of hydrodynamical 
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simulations of the formation of Wolf-Rayct shells 
like NGC 6888. The most complete of these is 
Garcia-Segura et al. (1996), based on a model for 
the evolution of a 35 M star. A summary of the 
stages in their scenario is as follows. 

During its main-sequence (MS) lifetime the star 
sweeps the surrounding interstellar material into a 
thin shell with a radius of tens of parsecs. The in- 
terior of the MS bubble is filled with hot (> 10 6 K) 
post-shock stellar wind material. In the absence of 
thermal evaporation of material at the outer edge, 
the electron density within the MS bubble should 
be of order 10~ 3 cm~ 3 . 

After 4.5 x 10 6 years the star evolves off the 
main sequence and into a red supcrgiant (RSG). 
Over its lifetime of 2 x 10 5 years , the RSG wind 
fills the inner few parsecs of the MS bubble with 
cold dense material. The characteristic mass loss 
rate of the RSG phase is ~ 10~ 4 M Q yr _1 with 
terminal velocity 15 km s _1 . Between the free- 
streaming RSG wind and the MS bubble a thin 
dense shell of post-shock RSG material forms. 

Lastly the star enters its Wolf-Rayet phase, 
with mass loss rate ~ 10~ 4,5 M Q yr -1 , with Voo ~ 
2000 km s^ 1 over its lifetime of 2 x 10 5 years. The 
WR wind progressively sweeps up the RSG ma- 
terial, eventually overtaking the outer RSG shell. 
The resulting collision fractures the outer shell and 
the WR bubble breaks out into the MS bubble 
material. At this point the nebular shell reaches 
maximum brightness due to its high density; ac- 
cording to their model, the nebula should be de- 
tectable over a period of order 10 years. 

The morphology and dynamics of WR shells 
are important not only as laboratories for study- 
ing wind interactions with the ISM, but also for 
the insights they provide into the mass-loss history 
of massive stars as well as the environments into 
which remnants of core collapse supernovae will 
propagate and interact. In this paper, we use ob- 
servations of a region of NGC 6888 obtained with 
the Hubble Space Telescope WFPC2 to study the 
ionization structure and infer physical conditions 
within the shell. 



2. Observations and Results 
2.1. Observations 

Observations of NGC 6888 were made in July 
1995 using the Wide Field/Planetary Camera 2 
on the Hubble Space Telescope. The orientation of 
the field of view is shown on a ground-based image 
of the nebula in Figure 1. The characteristics of 
the WFPC2 instrument are described in Trauger 
et al. (1994). Three of the cameras have an im- 
age scale of 0'.'0996 per pixel; the fourth camera 
has an image scale of 0'.'0455 per pixel. Two ex- 
posures were made in each of three filters: F502N 
([O III] A5007), F656N (Ha), and F673N ([S II] 
AA6717, 6731). The integration time for each fil- 
ter was 2 x 1100 seconds. Each pair of exposures 
was reduced and combined using a cosmic-ray re- 
jection routine. Remaining cosmetic defects were 
removed by direct inspection of the images. The 
resultant CCD frames were mosaicked using a rou- 
tine utilizing the astrometric solution in Holtz- 
mann et al. (1995a). Figure 2 shows a false- 
color composite of the narrow-band images in a 
linear stretch, with [O III] shown as blue, Ha as 
green, and [S II] as red. Grayscale images from 
the individual filters are shown in Figures 3 (Ha), 
4 QO III]), and 5 ([S II]). 

The conversion from DN to surface brightness 
(in erg cm -2 sec -1 pixel -1 ) for the narrow-band 
filters using the procedure outlined in Holtzmann 
et al. (1995b); 

F > = DNX (Md X A,Q E T xXtexp (1) 

where DN — total data numbers in a given pixel, 
GR t ~ 2 for the gain state used (7e"/DN), E x 
= energy of a photon at the observed wavelength, 
A = collecting area of the HST, QT\ = the tele- 
scope+filter efficiency at the observed wavelength, 
and t exp = the exposure time of the image. 

The flux must then be corrected for internal and 
interstellar extinction. Mitra (1990) estimates the 
extinction in the NW region from the H/3/Ha ratio 
in the spectra. Results for the three positions are 
Chp — 0.94, 1.05, and 1.13. Using the median 
value with the reddening law given by Cardelli et 
al. (1989), the dereddened surface brightness in 
this region is given by 
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/(F656A0 



and similarly, 



^(^656^) x 10 Cff ° 
5.27F(F6567V). (2) 



I(F673N) = 5.04F(F6737V), (3a) 
I(F502N) = 10.2F(F502iV). (3b) 



2.2. Results 

At a distance of 1450 pc one WF pixel corre- 
sponds to 2.18x 10 15 cm. The entire field of view is 
thus about 1 pc on a side. The network of clumps 
comprising the nebular shell are photoionized by 
the central star HD 192163, about 3 parsecs away 
in the direction indicated in Figure 2. We assume 
that the nebula is oblate; i.e., the clumps and cen- 
tral star lie approximately in the plane of the sky; 
hence, distances in the image do not suffer geo- 
metric foreshortening. Even accepting the prolate 
shape of Lozhinskaya (1970), distances would only 
be compressed by 15%. The photoionized clumps 
are visible to varying degrees in all three filters, 
with a strong correlation between [S II] brightness 
and Ha brightness. The filaments are typically 
shorter along the axis pointed towards the central 
star relative to the axis perpendicular to it. The 
network of filaments is enveloped by a "skin" of 
[O III] emission, formed by the shock driven into 
the external material by some combination of the 
nebular shell expansion and the pressure of the 
post-shock WR wind. The diffuse component of 
the nebular Ha and [S II] emission, mostly inter- 
nal to this skin, is a combination of emission from 
these two regions. In the NW quadrant, one clump 
deviates significantly from the positions of the rest 
of the shell. It has an [O III] arc extending beyond 
it, as well as a diffuse Ha tongue of emission ema- 
nating from the clump and extending beyond the 
arc itself. 

3. Photoionization Of The Nebular Shell 

The material lost by the central star in its pre- 
vious red supergiant wind has been swept up by 
the subsequent fast wind of the current WR phase. 



As the shell thinned it became subject to instabil- 
ity and fragmented. The bright patches of visible 
emission of NGC 6888 seen here are understood to 
be the fragmented circumstellar shell photoionized 
by the UV flux of the central star. 

3.1. Spatial Profiles 

Features were chosen on the basis of high Ha 
brightness. Those clumps with significant varia- 
tion at the resolution limit were rejected. Those 
clumps that appeared to have a complicated ge- 
ometry were also rejected. A total of six clumps 
across the observed region were chosen from the 
remaining sample. Spatial profiles were extracted 
along the direction towards the central star with a 
stepsize of 0.5 pixel. The positions of the spa- 
tial profiles are shown in Figure 6. The back- 
ground level for each profile was fit fit by a lin- 
ear or quadratic function and subtracted. Finally, 
the depth along the line of sight was estimated 
from the geometry of a given clump, using the ge- 
ometric mean of the major and minor axes for the 
spheroidal clumps and the overall curvature of the 
extended clumps as a rough equivalent to the emit- 
ting column. The volume emissivity of the gas is 
then given by 



4ttja = 



(4) 



where il = solid angle subtended by one WF pixel, 
and I = estimated line-of-sight depth, in cm. 

A first estimate of the hydrogen density can 
be derived from these cmissivities; using Oster- 
brock (1989), Table 4.2, one sees that for 10 4 K 

gas, 



Anjx = 3.56 x 10~ 25 



n p n e erg cm 3 s 1 



(5) 



For the clumps chosen the surface brightness 
and estimated line-of-sight depths result in densi- 
ties of order 1000 cm~ 3 . 

The spatial profiles extracted from the im- 
ages were modeled using the photoionization code 
CLOUDY (Ferland et al. 1998). Satisfactory fits 
were obtained by modeling clumps as thin slabs of 
constant density with two-component linear wings 
to the front and rear of the slab. The inflection 
point of the wings occurs where the density falls 
to one-half the density of the slab. Thus, each 
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clump model has five spatial parameters: the be- 
ginning and end of the slab, the two inflection 
points, and the total width of the clump model. 
The abundances in the gas were fixed at solar val- 
ues with the exception of helium, oxygen, nitro- 
gen, nitrogen, and sulfur. Helium was fixed at 
twice the solar abundance, following the results of 
Kwitter (1981). Varying this abundance did not 
significantly change the emissivity of the other el- 
ements. The other abundances were constrained 
by the observed profiles. For oxygen and sulfur, 
the elemental abundance was varied to maximize 
agreement with the appropriate observed profile. 
For nitrogen, the abundance was chosen to match 
the [N II] /Ha emission ratio to the value oberved 
by Parker (1978) for the corresponding area of 
the nebula. It is necessary to include the nitrogen 
emission in the model profile, as the F656N filter 
admits approximately the equivalent of 15% of the 
[N II]A6584 line flux (Birctta ct al. 1996), which 
for NGC 6888 can represent 20% of the counts in 
the F656N image. 

Each clump was placed at its projected distance 
from the central star. The input ionizing spectrum 
was the model atmosphere for a WN5/WN6 star 
presented in Hillier and Miller (1998). For com- 
parison to the observed profiles the model profiles 
were convolved with a gaussian of width corre- 
sponding to the PSF in the images, approximately 
0.5 pixel. 

3.2. Model Results 

A representative comparison of an observed and 
model profile is shown in Figure 7. Models with 
terminating ionization fronts consistently overes- 
timated the separation of the peak [S II] and Ha 
emissivities, suggesting that the shell does not con- 
tain large amounts of neutral material. This is 
supported by the absence of [O I] in the spec- 
trophotometry of Mitra (1990); based on observed 
line fluxes, I([0 I]A6300) < 0.02 I(H/3). The best 
fit photoionization models required a minumum H- 
ionizing flux of logQ = 49.2, with significant de- 
viations for two of the models (clumps 2 and 5) for 
values less than 49.3. This is in reasonable agree- 
ment with the value of 49.0 from the Wolf-Rayet 
models of Crowther and Smith (1996). It is, how- 
ever, at odds with Marston and Meaburn (1988), 
who calculate the Lyman continuum flux required 
to maintain the observed Ha brightness as logQo 



= 47.6. The implications of this discrepancy are 
addressed in §4. 

Models were constructed for logQo from 49.2 
to 49.4. The uncertainty on the abundance for an 
individual model is small; variation in the model 
elemental abundances by as little as 0.02 dex was 
sufficient to shift the emissivity of the correspond- 
ing line out of agreement with the observation. Us- 
ing the traditional nomenclature the abundances 
are (12 + log(N/H)) = 8.1 ±0.2, (12 + log(0/H)) 
= 7.8 ±0.2, and (12 + log(S/H)) = 6.8 ±0.2. The 
quoted uncertainty represents the spread between 
the abundances of all of the models, summarized 
in Table 1 and shown graphically in Figure 8. 

One of the primary sources of error in deter- 
mining abundances with this methodology is the 
uncertainty in the estimated linc-of-sight depth for 
a given clump. This somewhat counterintuitive re- 
sult arises from the fact that, while the ionic abun- 
dance can be inferred directly from the observa- 
tions independent of density or line-of-sight depth, 
the ionic fraction differs significantly with changes 
in the density and ionizing flux (i.e., changes in 
the ionization parameter). For example, a 25% er- 
ror in the estimated line-of-sight depth is roughly 
equivalent, in its effect on abundances, to an er- 
ror of about ±0.05 dex in Q . This typically cor- 
responds to an uncertainty of about ±0.1 dex in 
abundances. In addition, the offset between the 
peak [O III] and Ha emissivities were consistently 
underestimated by the models, which we ascribe 
to the simplicity of the one-dimensional profile 
model. This deviation suggests that our oxygen 
abundance may be systematically low; based on 
the model electron density at the position of the 
observed peak, this effect is at most 0.15 dex. 

3.3. Comparison with Ground-Based Spec- 
troscopy 

A longstanding question in nebular analysis in- 
volves the degree to which spectroscopy integrates 
over structure within nebulae, and the effect of 
inhomogeneity on inferences based on those data 
(e.g., the t 2 parameter in Peimbert 1967). Indeed, 
our images show significant density and ionization 
structure on spatial scales of 1" and shorter. A 
single resolution element typical of ground based 
spectroscopy would necessarily sample a range of 
physically diverse regions. Further, discrete fea- 
tures are typically embedded within diffuse line- 
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of-sight background that is comparably bright to 
the features themselves. It is unreasonable to 
imagine that ensemble-average conditions mea- 
sured spectroscopically have much physical mean- 
ing for specific locations within the nebula. Nor, 
conversely, can regions sampled spectroscopically 
be described meaningfully using a single set of 
physical parameters such as temperature or den- 
sity. 

An examination of inferred electron densities 
makes the point. Spectroscopy (Kwitter 1981; 
Mitra 1990) of this region typically yields [O II] 
and [S II] electron densities in the range from 100 
to 450 cm -3 , well below our densities of 1000-1600 
cm~ 3 . As discussed above, our density estimates 
are based on no more than the physics of hydro- 
gen recombination. Enforcing spectroscopically- 
detcrmined densities on features herein would re- 
quire that their line-of-sight depths exceed their 
projected sizes by one to two orders of magnitude. 
A simpler explanation is that ground-based spec- 
troscopy samples not only discrete features, but 
also unavoidably includes large volumes of lower 
density, more diffuse material. Roughly speaking, 
a mix of equal amounts of [S II] emission from gas 
at 1500 cm~ 3 (typical of our inferred densities) 
and [S II] emission from gas in the low density 
limit will yield a "measured" spectroscopic den- 
sity of around 500 cm -3 - even though there is 
little or no gas with that density actually present 
within the sampled volume. Comparable density 
inhomogeneities also exist within discrete features 
themselves. Similarly, our models show that T e 
can easily vary by 1000 K within the region sam- 
pled by a spectrograph aperture, even in the ab- 
sence of background emission. 

To employ traditional methods of interpreting 
spectra (e.g., icF$) it is necessary to assume a 
certain degree of homogeneity within a zone re- 
sponsible for emission from a given ion, as well 
as some correspondence between zones responsi- 
ble for emission from different ions. These ad hoc 
methods have a definite advantage - they yield an 
answer. However, the physical assumptions upon 
which they are based are not well supported by 
our data or our models. While we do not observe 
the wealth of lines present in a spectrum, we are 
able to use a few lines from both low and high ion- 
ization species to constrain realistic models of the 
ionization structure of observed features. Our re- 



sulting abundance determinations differ from the 
results of Kwitter (1981) by +0.3 dex, -0.7 dex, 
and and +0.5 dex for N, O, and S, respectively. 
The concerns that we raise are in addition to those 
of Alexander and Balick (1997), who found that 
the icf methodology may fail even in cases where 
the density is homogeneous. While the method 
we employ embodies its own suite of uncertain- 
ties, we suggest that there is ample evidence to 
warrant increased skepticism regarding the results 
of traditional nebular abundance determinations. 

4. Nature of the Skin Emission 

The dense shell is enveloped in a thin skin of 
emission, best seen at the outer edge of the nebula 
beyond the network of clumps, and most evident 
in [O III] . A perpendicular cut taken into the neb- 
ula at a bright part of the skin is shown in Figure 
9. Notice that there is a rise in [O III] emission, fol- 
lowed by a rise in Ha a few 10 16 cm closer towards 
the star. If the difference in ionization state across 
the shell were due to photoionization effects, we 
would expect that the higher ionization emission 
would be found interior to the lower ionization 
emission, as seen in the clump profiles discussed 
above. Rather, this morphology suggests instead 
that the leading edge of [O III] emission arises 
in the cooling behind a radiative shock driven 
into the relatively low density medium surround- 
ing the shell. The recombination region, how- 
ever, is "incomplete" , and the postshock flow cools 
and recombines only to the point that it comes 
into photoionization equilibrium with the ionizing 
flux from the star. Since the ionization parame- 
ter is high in this trailing photoionized zone, the 
[O III] /Ha ratio remains high (~ 5) and [S II] re- 
mains very weak throughout. There are, however, 
a few positions where post-shock material is more 
prominent in Ha. These are probably locations 
where the shock is partially "shadowed" from the 
stellar UV flux by dense clumps within the shell. 

An arc seen in the NW quadrant of the field 
in [O III] emission extends well beyond the gen- 
eral perimeter of the skin. This is morphologically 
consistent with the early stage of a "blowout" : a 
place where the shocked stellar wind has broken 
more completely through the nebular shell mate- 
rial and is expanding outward into the MS bubble 
itself. An echelle spectrum of this arc was taken in 
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the light of [O III] in December 1997 by Graham et 
al. (1998) using the 3.5-m telescope at Mt. Hamil- 
ton. We do not discuss this spectrum in detail here 
other than to note that it showed the character- 
istic velocity structure of an expanding bubble at 
this location. Assuming the line-of-sight angle of 
the expansion is the same as the angle of the arc 
to the general perimeter of the skin, the echclle 
arc has in an expansion velocity of 70 ± 12 km 
s -1 , implying a shock velocity of < 93 ± 16 km- 
s _1 . This is insufficient to produce such a large 
[O III] /Ha in a neutral pre-shock medium, sug- 
gesting that the pre-shock medium is fully ionized 
Cox and Raymond (1985). We show below that 
this is supported by observation. 

We estimate the emissivity of the post-shock 
gas at the leading edge of the NW arc, again using 
the geometry to estimate the pathlength. Assum- 
ing a temperature of 10 4 K, this yields a hydrogen 
density of 100 cm~ 3 . Due to Poisson noise, the 
pathlength estimate, and the unknown tempera- 
ture, the uncertainty is about a factor of 2. If we 
assume the pre-shock gas temperature is also 10 4 
K, then we can use the isothermal assumption to 
relate the densities of the pre-shock MS material 
and post-shock skin material; i.e., 

p skin = 2M 2 p M s- (6) 

With the shock velocity 93 km s" 1 , M = 5.6, and 
the resulting pre-shock density is 2 cm~ 3 . 

On the basis of our analysis, it is clear that 
the shell is extremely "leaky" to ionizing emis- 
sion. Using parameters derived below, the optical 
depth of the postshock region to ionizing radiation 
is only about 2 x 10~ 3 assuming logQo = 49.3, 
and so any radiation that reaches the shock has 
effectively escaped from the nebula. The densest 
knots themselves have optical depths of order 1. 
As can be seen from the more face-on portions of 
the shell in 1, such dense knots cover a very small 
fraction of the surface of the nebula. As a result, 
Zanstra methods cannot be used to obtain reliable 
estimates of the total H-ionizing flux of the cen- 
tral stars of objects such as NGC 6888. Based on 
the H-ionizing flux budget of the nebula obtained 
by Marston and Meaburn (1988) (logQ = 47.6) 
and the value derived in this work (49.3), only 2% 
of the ionizing photons from HD 192163 are pro- 
cessed within the nebular shell of NGC 6888. 



Another consequence of this analysis is that any 
material surrounding NGC 6888 will be fully ion- 
ized. This implies that there cannot be a signif- 
icant amount of neutral material close to the ob- 
served nebula as suggested by Marston (1995). 
However, an ionized shell of density 2 cm~ 3 ex- 
terior to the nebular shell could contain a signifi- 
cant amount of mass with low surface brightness. 
It may be that a major portion of the estimated 
18 M Q of material lost in the RSG wind has been 
lost to the MS bubble, and exists as a low density 
shell is only visible in post-shock emission. 

5. Discussion 

NGC 6888 is primarily composed of the mate- 
rial from the RSG wind of the Wolf-Rayet pro- 
genitor, swept into a shell by the central star's 
WR wind. This shell is photoionized by the UV 
flux from the central star. We combine the results 
of this paper and previous X-ray observations to 
create a coherent picture of the past and present 
state of NGC 6888. In Figure 10 is a sketch of the 
optical and X-ray emitting regions. 

5.1. Gas Pressure 

The gas pressure of a shocked stellar wind can 
be expressed in terms of its mass loss rate and 
its terminal velocity. For an adiabatic shock, the 
post-shock gas pressure is 



P - = Tex — E^T = 2 ' 88 x W 7 M- 4 V 3 R- 2 cm" 3 K, 

(7) 

where M_ 4 is the mass loss rate in units of 10~ 4 
M yr _1 , V3 is the terminal velocity in units of 10 3 
km s _1 , and R w is the radius of the termination 
shock in parsecs. 

Crowther and Smith (1996) derive the stellar 
parameters M_4 = 1.23 and V3 = 1.75 from model 
fits to infrared observations of WR 136. In addi- 
tion, they also summarize the results of previous 
studies for this object. If we assume that the inner 
edge of the X-ray emission corresponds to the lo- 
cation of the termination shock, we estimate that 
R w =3.1 pc. These results yield a pressure behind 
the stellar wind shock of 

P w = 3.4 - 5.3 x 10 6 cm~ 3 K. (8) 
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assuming that the location of the termination 
shock is not changing rapidly. 

The X-ray flux comes from the post-shock WR 
wind material. For the X-ray emitting region, the 
results of Wrigge et al. (1994) were used for the 
temperature, emission measure, and emitting vol- 
ume. (We note that we suspect the value quoted 
in their paper for the emission measure has a tran- 
scription error, resulting in an error of a factor of 
10 28 . We use this corrected value in this paper.) 
The density can be derived from the emission mea- 
sure E; i.e., 

n 2 H V ~ Ex AnD 2 , (9) 

where V is the volume of the emitting region, 
and D is the distance to the nebula. Wrigge et 
al. (1994) get an emitting volume of 7 - 12 x 10 55 
cm 3 and a (corrected) emission measure of (1.5 ± 
1.1) x 10 11 cm~ 5 , then at an assumed distance of 
1450 pc the density of the X-ray gas is n# = 0.3 - 
0.8 cm" 3 . With a temperature of 2.36 ±0.52 x 10 6 
K, the pressure of the X-ray gas is 

P x = 1.4-3.8 x 10 6 cuT 3 K, (10) 

roughly equal to that of the post-shock stellar 
wind. This gas is not simply the post-shock wind 
material, as its expected temperature is 4.2 x 10 7 
K and of such low density that it would be unable 
to cool quickly to the observed X-ray gas tem- 
perature. Instead, this material is probably the 
result of thermal evaporation of nebular shell ma- 
terial into the post-shock stellar wind (Hartquist 
et al. 1986). The rough equivalence between the 
pressure expected from eq. (8) and the pressure 
of the observed gas from cq. (10) supports this 
contention. 

The gas pressure of each of the regions can be 
calculated assuming an ideal gas law. For the 
outer shock and nebular shell, the values used are 
those derived in the previous sections. For the 
outer shock, the pressure is 

Pskin ~2x 10 6 cm~ 3 K. (11) 

Given uncertainties in parameters used to infer 
pressures for the various components of NGC 
6888, this is not significantly different from the 
pressures inferred for the shocked stellar wind and 
the X-ray gas. 



For the clumps comprising the nebular shell the 
central pressure can be calculated from the mod- 
els. Using a peak density of 1300 cm~ 3 at a tem- 
perature 10 4 K typical of the spatial profile models 
yields a pressure of 

Psheii = 26 x 10 6 cm~ 3 K. (12) 

Thus, clumps in the nebular shell are in severe 
overpressure compared to the surrounding gas, by 
a factor of order 10. This is consistent with the 
inclusion of density ramps in the model profiles of 
these regions, although the magnitude of the dis- 
crepancy is surprising. With typical clump size of 
10 17 cm, isothermal expansion of the clumps at 
the sound speed should alleviate this pressure dis- 
crepancy on a timescale of order 2000 years. This 
implies that until recently the termination shock 
was at Rm < 1.3 pc. If we assume that the loss of 
pressure support for the termination shock would 
result in its motion outward at the stellar wind 
velocity, it would take of order 1200 years for R^, 
to move from 1.2 to 3.1 pc. Considering that the 
clump pressure was initially higher than now ob- 
served coupled with the quadratic dependence of 
pressure with R^ , there can be a considerable time 
where the clumps have not yet come into pressure 
balance with its surroundings. 

Is there evidence supporting the contention of 
rapid depressurization of the bubble interior? Ini- 
tially, we note that the thickness of the postshock 
region behind the outer [O III] skin provides an es- 
timate of the age of this shock. Assuming a post- 
shock density of 100 cm -3 , a preshock density of 
2 cm~ 3 , a shock velocity of <~ 100 km s _1 , and 
a thickness of 10 16 cm, the shock broke through 
the densest portion of the shell about 1500 years 
ago. Blowouts can also be seen on the global scale 
of NGC 6888. In Miller and Chu (1993) Fig- 
ure 3b, one sees that the elliptical nebular shell 
is enclosed by a complete, approximately spheri- 
cal [O III] bubble. The loss of containment would 
cause a rapid drop of internal pressure like that im- 
plied by the pressure mismatch between the shell 
and stellar wind. If we use the time for the clumps 
to expand into balance with the new pressure as a 
constraint, then the edge of the blowout has to tra- 
verse 1.3 pc in much less than 2000 years, requiring 
an expansion velocity well over 600 km s . This 
implies that the density in the outer MS bubble 
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drops with increasing distance. Using the present 
internal pressure as a lower limit, the density of 
the outer MS bubble is still > 0.1 cm" 3 . 

5.2. Comparison with Hydrodynamical 
Simulations 

The most complete numerical simulations of 
Wolf-Rayet bubble formation and dynamics is 
Garcia-Segura et al. (1996). In it, they follow the 
evolutionary phases of a massive star, with main- 
sequence, red supergiant, and Wolf-Rayet mass 
loss. In summary, they produce an MS bubble 30 
pc in radius, filled with hot (10 7 K) thin (0.002 
cm~ 3 ) material; 18.5 M Q of RSG material, with a 
thin shell of post-shock RSG material at the RSG- 
MS interface; and a Wolf-Rayet wind that sweeps 
up the RSG material, leading to a collision with 
the thin RSG shell and eventually breakout of the 
WR wind into the MS bubble. The results of this 
paper show two main deviations from this picture. 
First, as argued before, the material filling the MS 
bubble is at least 0.1 cm~ 3 . For the MS bubble 
density implied by Garcia-Segura et al. (1996), 
any outer shock would not be visible. The thermal 
pressure of the observed material is the same as 
that assumed for the simulations, so a thin RSG 
shell can still form. Second, their post-collision 
shell density (1000 cm -3 ) is consistent with that 
derived in this paper, but the total mass of the 
shell is 4 times higher than that implied by our 
result of a fully ionized shell. A simple resolution 
to both of these differences is the presence of ther- 
mal evaporation from the thin RSG shell into the 
MS bubble. At 2 cm~ 3 , a ~ 1 pc thick annulus 
of material outside the nebular shell would have 
a mass of 15 M Q of material and still be unde- 
tectable in the WFPC2 image. A smaller annulus 
with a density gradient could hide a similar mass. 
It has also been suggested that the central star's 
proper motion can remove it from its relic MS 
bubble; this does not resolve the RSG mass loss 
issue. Another minor issue is that of the limitation 
of these simulations. In a non-spherical shell the 
pressure will drastically change before the shells 
collide along the major axis. How that alters the 
dynamics of the shell is undetermined, although 
the problem of ellipsoidal shells was pursued in 
Garcia-Segura and Mac Low (1995). 



6. Conclusions 

NGC 6888 is a bubble formed by the mass loss 
of the precursor to the central star WR 136. The 
internal pressure of the shocked WR stellar wind 
has swept up the ejecta from the previous RSG 
phase. This material collides with a thin RSG 
shell formed at the MS bubble interface, fragment- 
ing into the clumps seen today. Models of spatial 
profiles of selected clumps appear to be fully ion- 
ized, resulting in a log Qo > 49.2 for the central 
star, and average abundances for the shell material 
of [N, O, S] = [8.1, 7.8, 6.8]. The oxygen deple- 
tion appears even more severe than that seen in 
previous spectroscopic studies. 

The internal pressure of the stellar wind drives 
a shock into the material external to the visible 
shell. The shock is visible in [O III] and Ha as 
a skin enveloping the clump network. The den- 
sity of this material suggests that the MS bubble 
cooled efficiently. Combined with the mass of the 
visible shell compared to the total mass loss of the 
RSG phase, we argue that thermal evaporation 
of material from the outer RSG shell has already 
occurred. Comparing the pressure inside and out- 
side the nebular shell suggests that a radical drop 
in pressure has recently occurred, and that the 
clumps have not yet come into equilibrium with 
the current pressure of the stellar wind. 
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Fig. 1— Ground-based image of NGC 6888, 
showing the field observed with WFPC2. North 
is at the top and east is to the left. The WFPC2 
footprint is approximately 150" (1 pc at a distance 
of 1450 pc) on a side. 



Fig. 2. — False-color composite of the WFPC2 im- 
ages of NGC 6888. Emission from [O III] is shown 
in blue. Emission from Ha is shown in green. 
Emission from [S II] is shown in red. 

Fig. 3.— WFPC2 image of NE section of NGC 
6888 taken through the F656N filter, which iso- 
lates light from Ha A6563. North and east are 
indicated by the arrows. The bar has a length of 
10", corresponding to 2.2 x 10 17 cm. The negative 
display uses a square-root stretch to allow com- 
parison over a large range in surface brightness. 



Fig. 9. — Spatial profile of the outer shock of NGC 
6888. The Ha surface brightness has been in- 
creased by a factor of 5 to ease comparison; [S II] 
has been raised by a factor of 10. Distance in- 
creases toward the central star. 



Fig. 4.— WFPC2 image of NE section of NGC 
6888 taken through the F502N filter, which iso- 
lates light from [O III] A5007. The negative dis- 
play uses a square-root stretch. 



Fig. 5.— WFPC2 image of NE section of NGC 
6888 taken through the F673N filter, which iso- 
lates light from [S II] AA6717,6731. The negative 
display uses a square-root stretch. 

Fig. 6. — Positions and widths of spatial profiles 
extracted from the WFPC2 data, superimposed 
on the Ha image. 

Fig. 7. — a. Volume emissivities from spatial pro- 
file 3. The distance into the clump is measured 
from the point towards the star at which the Ha 
profile flattened out. Ha is plotted with a solid 
line; S II, dotted line; O III, dashed line. b. Pho- 
toionization model for log Qo = 49.3 that best 
matched the observed profile, c. The data (solid) 
with the model (dotted) overplottcd. Distance in- 
creases away from the central star. 



Fig. 10. — a. Sketch of the wind-blown bubble as 
the red supcrgiant (RSG) wind is being swept up 
by the Wolf-Rayet (WR) wind into a "WR shell" . 
The "MS bubble" is composed of a mixture of 
relic main-sequence wind material and evaporated 
"RSG shell" material, b. The bubble some time 
after the WR and RSG shells collide. The post- 
shock WR wind now streams past the fragmented 
RSG shell and directly interacts with the MS bub- 
ble, driving a shock visible as a "skin" enveloping 
the nebular shell. The MS bubble material is shad- 
owed from the UV flux of the central star by the 
densest of the knots in the nebular shell, affecting 
its [O III] flux. 



Fig. 8. — Plot of model nitrogen, oxygen, and 
sulfur abundances as a function of the total H- 
ionizing flux log Qo- The lines are identified by 
their correspondent spatial profile number. 
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Table 1 

Model Parameters for Spatial Profile Fits 



log Qo 



Spatial Profile 






49.20 


49.25 


49.30 


49.35 


49.40 


1 


Peak iih (cm 


) 


1570 


1570 


1540 


1530 


1520 




N 




8.02 


8.02 


8.06 


8.11 


8.18 









7.77 


7.71 


7.67 


7.63 


7.59 




s 




6.45 


6.52 


6.60 


6.68 


6.71 


2 


Peak nH(cm" 




1300 


1310 


1320 


1320 


1320 




N 




8.03 


8.04 


8.05 


8.09 


8.16 









7.68 


7.57 


7.57 


7.52 


7.47 




s 




6.69 


6.79 


6.81 


6.87 


6.91 


3 


Peak iih (cm~ 


- 3 ) 


1030 


1020 


1010 


1000 


1000 




N 




8.11 


8.18 


8.26 


8.34 


8.43 









7.54 


7.48 


7.47 


7.41 


7.40 




s 




6.72 


6.82 


6.90 


6.92 


6.98 


4 


Peak iih (crrT 


- 3 ) 


1150 


1140 


1140 


1140 


1130 




N 




7.69 


7.72 


7.77 


7.85 


7.95 









8.05 


7.98 


7.89 


7.86 


7.82 




s 




6.68 


6.77 


6.83 


6.89 


6.96 


5 


Peak iih (cm" 


- 3 ) 


1100 


1100 


1100 


1130 


1120 




N 




8.13 


8.15 


8.16 


8.19 


8.26 









7.77 


7.77 


7.72 


7.72 


7.67 




s 




6.53 


6.56 


6.61 


6.66 


6.76 


6 


Peak iih (crrT 




1500 


1500 


1490 


1490 


1480 




N 




7.96 


7.98 


8.03 


8.10 


8.19 









7.87 


7.77 


7.72 


7.65 


7.61 




s 




6.89 


6.94 


7.00 


7.04 


7.08 
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